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tJ" . Abstract 
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I We investigate the dependence of hadronic ratios on the number 

■ of projectile participants using a thermal model incorporating exact 

. baryon number and strangeness conservation. A comparison is made 

with results from Au — Au collisions obtained at the BNL-AGS. 
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Prehminary results on the dependence of hadronic ratios on the number 
^ \ of projectile participants have recently been presented by the E866 collabo- 

^ ' ration |Q for relativistic Au — Au collisions at the BNL-AGS. These results 

give insight into the behaviour of the produced hadronic system as a function 
^ ■ of the baryon number and of the size of the interaction volume. 

It is the purpose of the present paper to analyze these results using a 
thermal resonance gas model at a fixed temperature and a fixed baryon den- 
sity. Our treatment differs from previous 0, |^ ones in that we consider the 
baryon content exactly. This means that we do not introduce chemical po- 
tentials for the baryon number (nor for strangeness). Chemical potentials 
are usually introduced to enforce the right quantum numbers of the system 
in an average sense. This is a correct treatment for a large system, however, 
for a small system the production of e.g. an extra proton - anti-proton pair 
will clearly be more suppressed than in a large system. These extra correc- 
tions were first pointed out by Hagedorn |^ and subsequently a complete 
treatment was presented by many people 0, H, |^, ||. We emphasize that 
these corrections do not contain any information about the dynamics. They 
simply follow from baryon number conservation. These corrections must be 



taken into account before considering more involved models. It is also worth 
emphasizing that they do not involve any new parameters. 

As an example we analyze the (preliminary) data recently presented by 
the E866 collaboration at BNL. 

The exact treatment of quantum numbers in statistical mechanics is ob- 
tained by projecting the partition function onto the desired values of B and 
S 

Zb,s = r e-^^<^ ^ r di,e-'^^Z{T, Xb, Xs) (1) 
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where the usual fugacity factors Xb and Xs have been replaced by : 

Xb = e''^ Xs = e'^. (2) 

We will use 

B = 2Npp (3) 

where Npp is the number of projectile participants with the factor 2 reflect- 
ing the symmetry of the Au — Au collision system. As the contributions 
always come pairwise for particle and anti-particle the fugacity factors will 
give rise to the cosine of the angle. In the further treatment it is useful 



to group all particles appearing in the Particle Data Booklet |T0[ into four 
categories depending on their quantum numbers (we leave out charm and 
bottom). Zk is the sum (given below) of all mesons having strangeness ±1 
{K, K, K*, . . .), similarly is the sum of all baryons and anti-baryons hav- 
ing zero strangeness, Zy is the sum of all hyperons and anti-hyperons while 
Zq is the sum of all non-strange mesons, and so on : 

d^p _Ej/T 



je\s\=i,\B\=o 



ie|5|=o,|B|=i 



ie|5|=i,|B|=i 



(27r) = 



^0 = E yf^s'-''^^- (4) 

je\S\=0,\B\=0 ^ W 

We do not include cascade particles as their contribution is unimportant 
for the energy range under consideration and their inclusion considerably 
complicates the formalism. Each term will be multiplied by the cosine of an 
angle, either or ip, in the case where two angles are needed (e.g. for the 
hyperons) one introduces a new one, a, using 

oo 

5{<p-i)-a)= e*"('^-^-"). (5) 
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Using the integral representation of the modified Bessel functions 



4(2) = - / e''°''' cos n9 d9 
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one can write the partition function as 

Zb,s = Zq In+B 



(6) 
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In order to discuss the particle abundances it is useful to introduce the fol- 
lowing quantities M 
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E ^n{2Zy)In+Bi2Zi\f)In+l{2ZK), 
n=—oo 
00 

E ^n+li2Zy)In+B~l{2ZM)In{2ZK), 

= — OD 
00 

E -^n+l(2^y)-^n+_B(22'7v)/„(2Zi<'), 

=—00 
00 

E In{2Zy)In+B{2Zi\f)In-l{2ZK), 

=—00 
00 

E -^n+l(2^y)/n+B+l(2^Ar)/n(2^_ft:), 
=—00 
00 

E -^n-l(2^y)/n+B(2^Ar)/„(2Zi^). 
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If a particle, i, has strangeness 1 and baryon number 0, it's density will 
be given by 



Zo 



R 



K 



'B,S 



(2^ 



-E,/T 
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while a particle with strangeness and baryon number 1, will have a density 
given by 



R 



N 
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-Ei/T 



(10) 



All other particle densities are obtained by using the appropriate R factor 
given in equation (8). The factor in square brackets replaces the fugacity in 
the usual grand canonical ensemble treatment ^ . Having thus determined 
all particle densities, we consider the behaviour at freeze-out time. In this 
case all the resonances in the gas are allowed to decay into lighter stable 
particles. This means that each particle density is multiplied with its appro- 
priate branching ratio (indicated by Br below). The abundances of particles 
in the final state are thus determined by : 



nK+ = ^riiBrii K^), 
n^- = ^niBr{i tt"), 
riK- =^niBr{i K~), 
np = YjiT-iBrii p), 

Up = ^niBr{i ^ p). (11) 

where each sum runs over all particles contained in the hadronic gas. 

The comparison with experimental results is shown in figures 1 to 4. To 
compare with earlier calculations [0, ^ we keep the temperature T and the 
baryon density B/V fixed. This corresponds to keeping the baryon chemical 
potential fixed in the standard hadronic gas calculations using the grand 
canonical ensemble. 

In figure 1 we compare our results with recent data from the AGS ^ |TT| . 
Figure 1 shows the /tt^ ratio. As one can see the results obtained from 
our calculation show a steep rise with Npp before leveling off. The dependence 
on the baryon density is minimal in this case. This result is confirmed by 
calculations done in the grand canonical ensemble which also show that this 
ratio is almost independent of the baryon density 0. We note that the 
experimental data indicate a slower rise than the model calculation. 

In figures 2, 3 and 4 we show the K~ /tt^ , the / K~ and the p/yr^ 
ratios. In each case good agreement is obtained with the results of the E866 
collaboration |T]|. The relevant temperature is around T k, lOOMeV, the 
baryon density is in the range B/V ~ 0.02 - 0.05 fm~^. In the grand canon- 
ical ensemble this corresponds to a baryon chemical potential of ~ 540 
MeV. 
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Figure Captions. 

Figure 1 : The /t:~^ ratio as a function of the number of projectile par- 
ticipants Npp. The sohd hue is obtained for T = 96 MeV and B/V = 0.024 
fm~^, the dashed hue corresponds to T = 103 MeV and B/V = 0.050 fm~^ 
while the dotted line corresponds to T = 100 MeV and B/V = 0.04 fm-=*. 
Figure 2 : The K~ /t^^ ratio as a function of the number of projectile par- 
ticipants Npp. The notation is the same as in figure 1. 

Figure 3 : The / ratio as a function of the number of projectile par- 
ticipants Npp. The notation is the same as in figure 1. 

Figure 4 : The p/tt"*" ratio as a function of the number of projectile partici- 
pants Npp. The notation is the same as in figure 1. 
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